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Abstract 



Existence of minimal length is suggested in any quantum theory of gravity such as string theory, 
double special relativity and black hole physics. One way to impose minimal length is deforming 
Heisenberg algebra in phase space which is called Generalized Uncertainty Principle (GUP). In 
this paper, we develop statistical mechanics in GUP framework. Our method is quite general and 
. does not need to fix the generalized coordinates and momenta. We define general transformation 

in phase space which transforms usual Heisenberg algebra to a deformed one. In this method, 
quantum gravity effects only acts on the structure of phase space and we relate these effects to the 
density of states. We find an interesting phenomenon in Maxwell-Boltzmann statistics which has 
not a classical analogy. We show that there is an upper bound for the number of excited particles 
in the limit of high temperature which implies to the condensation. Also we study modification of 
\ Bose-Einstein condensation and the completely degenerate gas. 



PACS numbers: 04.60.Bc, 05.20.-y, 05.30.-d 
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1 Introduction 



In Hamiltonian viewpoint of classical mechanics, canonical equations of motion can be represented by 
the Poisson brackets in which the coordinates X{ and their conjugate momenta pj obey the Poisson 
algebra {x{,pj} = 5ij. Physically, these relations may be interpreted as the measurement of the 
position and momentum of a moving particle simultaneously. This means that the coordinates of each 
point of the corresponding phase space can be determined without any uncertainty. The transition 
to the quantum mechanics is straightforward. The classical dynamical variables should be replaced 
by their hermitian operators counterparts in Hilbert space and the Poisson brackets with the Dirac 
commutators. Hence, the above Poisson algebra takes the form of the Heisenberg algebra = 
ihSij. A consequence of such a algebra is that it divides the phase space to pixels proportional to 
the Planck constant and thus, the position and the momentum of particles can not be determined 
simultaneously. This is what we know as the uncertainty principle in the usual quantum mechanics. 
On the other hand, one of the most important predictions of the theories which deal with quantum 
gravity is that there exists a minimal length below which no other length can be observed [1]. One 
of the interesting features of the existence of such a minimum length is the modification it makes to 
the standard commutation relation between position and momentum in ordinary quantum mechanics 
which is called the generalized uncertainty principle (GUP) [2]. In one dimension the simplest form 
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of such relations can be written as 

AAAP>^(l + /3(AP) 2 + 7 ), (1) 

where j3 and 7 are positive and independent of AX and AP, but may in general depend on the 
expectation values < X > and < P >. The usual Heisenberg uncertainty relation can be recovered 
in the limit /3 = 7 = 0. It is easy to see that this equation implies a minimum position uncertainty of 
(AX) m i n = Ti\f]3. For a more general discussion on such deformed Heisenberg algebras, especially in 
three dimensions, see [3]. Now, it is possible to realize equation (1) from the following commutation 
relation between position and momentum operators 

'Xi,p 3 ] = in{i + pp 2 )Sij, (2) 

where P 2 = PiP 1 and we take 7 = (3 < P > 2 . Also, assuming that 

'Pi,Pj]=0, (3) 
the commutation relations for the coordinates are obtained as 



Xi , Xj 



= 2iHP (PiXj - PjXi), (4) 



which means that in more than one dimension, GUP naturally implies a noncommutative geometric 
generalization of position space. 

In a statistical mechanics point of view, the microstate of a given classical system may be defined 
by 3N position coordinates X±, X^n and SN momenta coordinates P±, P^n, where iV is the 
number of particles in the system. In a geometric picture, the set of coordinates pQ,P,), where 
i = l, ...,3iV, may be considered as a point in a 6N dimensional space, the so-called phase space of 
the system. Since the coordinates Xi and Pi vary with time, the dynamics of the whole system can be 
determined with the help of the Hamiltonian equations of motion for each of these coordinates, that 
is, Xi = {Xi,H} and Pi = {Pi,H}, where H(Xi,Pj) is the Hamiltonian of the system. Therefore, 
the dynamical behavior of the system may be viewed as a continuous trajectory of the phase point 
(Xi,Pi) in the phase space. On the other hand, it is well known that the correct description of 
such a system and its evolution requires quantum mechanics. In this case, the quantum particle 
has no well-defined trajectory in the phase space. This is because that at a given time t, there is 
no longer an exact position, but a typical extent of the wave function inside which this particle can 
be found. Another way of understanding this issue is to state the Heisenberg uncertainty principle 
which is a result of the commutation relations between position and momentum operators. In this 
sense, if the position coordinate Xi is known up to AXj, then its conjugate momentum Pj cannot 
be determined with an accuracy better than APj, such that AAj A Pj > fi/2. The phase space is 
therefore divided into cells of area of the order of h 3N . The accuracy on the phase space trajectory is 
thus limited and this produces a discretization of the phase space. Now, what happened if one takes 
into account the GUP considerations instead of the ordinary uncertainty principle? The motivation 
of this question is that the GUP scheme relies on a modification of the canonical prescriptions and, 
in this respect, can be reliably applied to any dynamical system. Over the past few years, a number 
of works have been done, in the area of statistical mechanics in the GUP framework. For instance, 
the thermodynamics of the ideal gas and ultra-relativistic gas in micro-canonical ensemble in the 
GUP framework are studied in [4]. For harmonic oscillators and ideal gases in canonical ensembles 
see [5]. The deformed density matrix is studied in [6] and modified uncertainty relation for inverse 
temperature and internal energy is addressed in [7]. Black body radiation with minimal length effects 
is considered in [8]. For thermodynamics in doubly special relativity theory which also predicts a 
minimal length, see [9]. Thermodynamics in the framework of non-commutative Heisenberg algebra 
(q-deformed) is investigated in [10]. 
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In this paper, our aim is to study the aspects related to the application of GUP in the framework 
of classical and quantum statistical mechanics. Our study covers the thermodynamics of the fermionic 
and bosonic systems when their underlying canonical structures are of the type of GUP-deformed. 
We deal with the density of deformed phase space and study the effects of existence of minimal length 
on the density of states. We show that there exist a general transformation which transforms the 
usual Heisenberg algebra to a deformed one and after we convince ourselves to the existence of such a 
transformation, we obtain a density of state with the help of the Jacobian of the transformation. Then, 
we study the thermodynamics of the ideal gas and extreme relativistic gas with the help of density of 
deformed phase space. The Bose-Einstein condensation and the completely degenerate gas are the case 
studies which we will deal with them in the GUP framework. Also, we find an interesting phenomenon 
in the Maxwell-Boltzmann statistics which we interpret it as the gravitational condensation and we 
show that it occurs for the ideal gas and the extreme relativistic gas in the same manner. 

2 Statistical mechanics 

Usually introducing the GUP version of a quantum theory is done by replacing the ordinary com- 
mutation relations by their generalized counterparts (2) and (4). However, in what follows we act 
differently and define a transformation in the phase space which transforms the usual Heisenberg 
algebra to a deformed one. In this picture, GUP effects only change the density of the phase space 
of a statistical system. As we will see the volume of phase space becomes larger while its degrees of 
freedom decreases. 

The microstates of any physical system are determined with quantum mechanics and its energy 
levels should be obtained from the Schrodinger equation. In GUP framework, the Schrodinger equa- 
tion becomes a non- linear or higher order differential equation and it is not easy to solve it. For 
example, for the wave function and energy spectrum of harmonic oscillator see [2] and [11]. Partic- 
ular non-linear Schrodinger equation in GUP framework is proposed in [12] and solved for quantum 
bouncer in [13]. For the higher order modified schrodinger equation for quantum mechanical systems, 
see [14]. Of course, solving the higher order differential equation has some dramatic consequences 
because of the mathematical difficulty due to non-linearity and initial values for higher order dif- 
ferential equations. In our method it is not necessary to solve the modified Schrodinger equation, 
instead, we work in semi-classical approximation and we use the classical Hamiltonian in deformed 
phase space. In quantum picture, energy of a given macrostate is a summation of the energies of all 
corresponding microstates. On the other hand, energy of the macrostate is given by the Hamiltonian 
which is a continuous function of the phase space variables. It is customary to approximate the 
summation over the energies of the microstates by the integral over all phase space variables. As we 
have mentioned, the phase space is constructed by the fundamental cells as order of Planck constant 
h, [15]. We show that the fundamental cell becomes larger and also has a momentum dependence 
as h GUP = h (1 + (3 P 2 ). Then, we can approximate summation over energies of microstates by 
integral over all deformed phase space with spotting lattice structure of phase space. Let us now 
explain our approach to transition from ordinary phase space to deformed phase space and also show 
the momentum dependence of the h GUP . In ordinary quantum mechanics, the coordinates X{ and 
momenta pi satisfy Heisenberg algebra [xi , pj] = ihSij where its classical counterpart is the Poisson 
algebra {x^ , pj} = . In the GUP framework, these relations take the form 



{Xi,Xj} 



2(3{P l X J - PjXi) 



(5) 



{Xi,Pj} 



(1 + /3P 2 )<%, 



(6) 



{Pi,Pj} =o. 



(7) 
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According to the Darboux theorem [16], it is always possible to find canonically conjugate variables 
Xi(x,p) and Pi(x,p) such that they satisfy relation (6) and (7). However, we would like to find 
variables Xi and Pi which satisfy also relation (5) simultaneously. In Appendix A, we convince 
ourselves that it is always possible to find generalized variables Xi and Pi which satisfy relations (5), 
(6) and (7) simultaneously. This possibility allows us to define general transformation 



(Xi,pi) ->• (Xi(x,p), Pi{x,p)), (8) 

which transforms the usual Heisenberg algebra to a deformed one. Now, we only should have the 
Jacobian of this transformation for our study. The Jacobian J of the transformation can be expressed 
in terms of poisson brackets in 2AA-dimensional phase space ( [5], see also Appendix A), 

3 = 4?^r = IK* ■ ^ = 0- + ?P 2 f- (9) 

d(xi,Pi) ^ 

In Appendix A we show that the relation (5) does not contribute in Jacobian, thus we miss the effects 
of non-commutativity in our study. Having the Jacobian of the transformation at hand, we can define 
the state density a(e) de in GUP framework by which we may approximate the summation of energy 
spectrums with an integral over deformed phase space. In the limit of large volume V, we can replace 
summation over energy spectrums with an integral as ~ * I zf~-> where dw = cf x d^p and M is the 
number of degrees of freedom. Here, ujq is a fundamental volume and is given by ujq = ', where h is 
the Planck constant. We define ordinary density of state as J2 — * I ao(e)de. Now, we can define the 
density of states in GUP framework with the help of transformation (8). In terms of the generalized 
variables Xi and Pi, the density of states for the deformed phase space becomes 



I cfXcfP cfXcfP 

a{£) d£ = w—j— = m^ww • (10) 

From the above relation, it is clear that the fundamental volume element of the corresponding phase 
space takes the form h(l + f3P 2 ) which is larger in comparison with the ordinary case. Thus, the 
number of accessible microstates for the system will decrease in GUP framework. Relation (10) plays 
an important role in our approach and we formulate the ensemble theory in GUP framework with 
the help of it in this section. 



2.1 Micro-canonical ensemble 

In equilibrium statistical mechanics, a system is supposed to be in the thermal equilibrium and 
different ensembles have different density of microstates in phase space. In pervious section we 
obtain density of the deformed phase space, then it is easy to generalize ensemble density for any 
ensemble in GUP framework. In micro-canonical, the energy of a given macrostate remains constant 
along the thermodynamical processes. This condition imposes the energy constraint to the system. 
Under this condition the probability of all microstates is same and density of state is constant. We 
assume that the density of states to be p(x,p) mce oc S(H(x,p) — E) (see [17]), where 5 is the Dirac 
delta function, which clearly satisfies two conditions, the same probability for all microstates and 
energy constraint. Now we define the number of microstates for N particles in GUP framework as 

j3N V ,3N _ N 

where E 1 = E'((3, E) is energy of each particle in GUP framework 1 . It is clear that the number 
of accessible microstates decrease because of the GUP correction in denominator which may be 

1 We work in unit where k b = 1 = c , where k b is the Boltzmann constant and c is the speed of light. 
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interpreted as a lost of information (see [18] and [4]). Now, the entropy of the system can be easily 
obtained from its standard definition as 



S=ln£(/3). (12) 

There is an interesting case, because of the decreasing number of microstates, entropy decreases and 
tends to zero in high energy regime (high temperature limit). In the very early theory of universe 
zero entropy coincides with preliminary singularity of universe (highest level of energy) and all the 
information miss there as we expect from the big-bang theory. Thermodynamics of micro-canonical 
ensemble can be obtained from entropy (12). 

2.2 Canonical ensemble 

In canonical ensemble the energy of a system is variable and system can exchange energy. It means 
that the energy of the microstates determine the probability of them in phase space. Then, the 
density of microstates should be a function of energy and we know it is proportional to Boltzmann 
factor: p(x,p) ce oc exp ( — H ^' p ^ ) . The partition function of the system is given by 

Z=£exp(~), ( i3) 

where e is the energy of the microstates and is a solution of Schrodinger equation. In fact we should 
solve modified Schrodinger equation in GUP framework and substitute energy of microstates in (13). 
However, in our approach we approximate summation over energy levels by an integral with the help 
of relation (10). Therefore, for a system with N particles we get 

rr , m ,r ^ 1 f f ct N Xd iN P . H(X,P). 

*„CT, V.ffl = w j J {1+pp2r exp( - -LM). (14) 

Now, all of the thermodynamical quantities in GUP framework can be obtain from the partition 
function [5]. 

2.3 Grand canonical ensemble 

In grand canonical ensemble, in addition of its energy, the system can be exchange its particles. Then 
we expect the probability of microstates depend to the number of particles and energy. As same as 
the pervious section the grand partition function in GUP framework is 



E(T,V,N,P)=Y,z N Z N (T,V,p), (15) 

JV=0 

where z = e M//T is the fugacity of the system and /i is the chemical potential. Again all the thermo- 
dynamical quantities can be obtained from grand canonical partition function. On the other hand, 
one may work with mean occupation number, < n £ > instead of grand partition function (15) which 
can be defined as 

< n e > = - -= , (16) 

z e 6 ' 1 + e 

where e is a constant with the values e = 1 (for Fermi-Dirac statistics) and e = — 1 (for Bose-Einstein 
statistics). The particle density and pressure can be obtain from < n £ > by definitions 

N = £ < n £ > , (17) 
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^ = ;El*(l + -^ /T )- (18) 

In the limit of Large volume, we can replace the above summation by an integral over all deformed 
phase space with the help of density of state in GUP framework (10). Fortunately, there is a closed 
form for the solution in Maxwell-Boltzmann statistics correspond to the limit z — > 0. In this limit we 
have a classical limit of minimal length and we only neglect the quantum effects, but the GUP effects 
still exist. On the other hand, in fully quantum area (z — > 1), we expect to find the purely quantum 
phenomenon such as Bose-Einstein condensation. In these limit, there is an analytical solution for 
any value of z in term of Fermi-Dirac and Bose-Einstein functions (see Appendix B). We study some 
consequences of minimal length for two gaseous systems, Ideal gas and Extreme relativistic gas in next 
section and consider thermodynamics of them in Maxwell-Boltzmann, Fermi-Dirac and Bose-Einstein 
statistics. 

3 Maxwell-Boltzmann statistics 

To study the effects of minimal length in thermodynamics, we can exert GUP condition with the help 
of relation (10), but we should note when we approximate summation in relation (17) by integral, 
the contribution of the first term which corresponds to the ground state energy of particles (e = 0) 
missed. Thus, we decompose the first term as 



N e = N — Nq = J <n £ > a(e) de, (19) 

where N e and Nq denote the number of particles in the excited states ground state respectively. The 
value of iVo which corresponds to y£- is negligible for classical limit z — > and becomes important 
in the limit of fully quantum area z — > 1. However, we keep this term in both cases. In Maxwell- 
Boltzmann statistics, the quantum effects are negligible and the type of particles (fermion or boson) 
is not important. Then, our start point is relations (17) and (18) according to which we have 



N d = Um ^ < n * > ) = z S e ~ £/T ' ( 2 °) 

and 

^ = Hm Q £ ln(l + eze-^ = z £ e~^\ (21) 

which resemble Maxwell-Boltzmann distribution as we expected for classical systems. It is seen that 
both of the above relations have the same value in this limit. In fact, they are a result of equation of 
state in classical statistical mechanics and by combining them, we get 



■pcly = N cl T ^ ^ 

which is a familiar classical equation of state. Thus the form of equation of states preserve for 
statistical systems in GUP theory. It is important to realize that the pressure V cl and N cl are in 
GUP framework and have different forms in comparison with their ordinary cases. Now, we calculate 
V cl and N cl for two gaseous systems, the ideal gas and extreme relativistic gas. 
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3.1 Ideal gas 



We consider a gaseous system consisting of N non-interacting particles confined in volume V which 
can sharing energy together. The energy momentum relation for ideal gas is given by e = where 
m is the mass of each particle. Then for sufficient large volume the summation in relation (20) can 
be replaced by an integral over all phase space with the help of density of states in GUP framework 
(10). By substituting (10) in (20), we have 



»«•/ f r d^Xd 3 P , 

N cl ~ z / / e ~2^T (23) 

6 " Jy Jo HI + PP 2 )f ( ' 

The integral over coordinates gives the three dimensional volume element V and we get 



N 



cl 



VZ 1^3/2 

4/i 3 V 



7 7 V a/7 



(24) 



where 7 = f3mT and erf is the error function. Here Nf denotes the number of excited particles for 
ideal gas in gravitational statistics. We can remove the gravity effects in the limit j3 — > 



lim Nf = ^z, (25) 
/3->-o A 3 

where A = , h is the mean thermal wavelength, in agreement with usual statistical physics. 
Relation (24) has a interesting consequence in the limit of T — > 00 as 



lim Nf = 

T-s.00 



VZ 7T 2 
4/l 3 03/2 ' 



(26) 



The above relation shows that there is an upper bound for the number of excited particles in the 
limit of very high temperature which is a purely quantum gravity effect and does not have a classical 
analogy. More precise, chemical potential is the mean energy which one particle need to apply to 
statistical system, then if we increase the number of particles, the energy of the system becomes 
large and larger. In ordinary Maxwell-Boltzmann statistics, the particles access to the upper energy 
levels in the limit of high temperature and there is not any limitation for them. However in GUP 
theory, there is a highest energy level for systems. Then, after the energy of the system reaches to 
the highest level, if we increase the number of particles, these should be lie in ground state and we 
have a condensation for classical ideal gas. Such a condensation is different from ordinary ones in 
Maxwell-Boltzmann statistics which is occurs in the limit of low temperature. By having Nf from 
relation (24), the pressure of the system can be obtained from relation (22) as 



~>cl 



Z IT 



3/2 1 



4mh 3 /3 5 /2 



( 7 2 -2 7 -l) 



7 



7 



erf(— ) 

V7 



(27) 



Also, the energy of the system can be obtained by definition 



U = T 2 

V 



-(-) 



(28) 



from which we get 



7 



u cl = 



zVtt 3 / 2 1 
8m/i 3 /?V2 7 2 



20Fe 1/7 ^/7 ( 7 2 + 4 7 + l)erfc( 



-7 ( (V2 - 2) 7 2 + (3^2 + 4) 7 + 2 ) 



(29) 



where erfc is the complementary error function. The pressure and energy of the classical ideal gas 
in purely classical limit (quantum effects are negligible) do not show the unusual behavior and have 
ordinary interpretation. 

3.2 Extreme relativistic gas 

As same as pervious section we consider N non-interacting particles which confined in volume V, but 
here the energy momentum relation for extreme relativistic gas is given by e = P . For large volume, 
we can replace summation in (20) with the help of density of state (10) and we have 



cl = Vz_ VtF 3 1 

h3 (]3/2 ^00 



4/3T 2 




(30) 



where G is the Meijer function. Like the ideal gas, we find a condensation due to the GUP theory 
for extreme relativistic gas in the limit of high temperature 



V 7 7T 2 

lim N? = 773-^. (31) 

Thus, again there is an upper limit for the number of particles in the limit of high temperature 
T — > 00. The maximum value of the number of particles for extreme relativistic gas is as same as the 
ideal gas (see relation (26)). The pressure can be obtained by substituting (30) in (22). For energy 
U cl of the system we have 



U 



cl 



IT ZV 3 1 



4/3T 2 



\ 



(32) 



In the limit of high temperature T — > 00, for the above relation we have 



U d - — — (33) 

Thus, there is an upper bound for energy of the extreme relativistic gases in GUP framework which 

is as order of Planck energy U cl ~ E pl (note /3 ~ in our unit). 

P i 



4 Fermi-Dirac and Bose-Einstein statistics 

In this limit the type of particles, which may be fermion or boson, plays an important role, according 
to which there are two types of statistics the so-called Fermi-Dirac and Bose-Einstein for fermions 
and bosons respectively. To consider the thermodynamics of fermions and bosons, we start again 
with relations (17) and (18). We replace summation in these relations by integral over all deformed 
phase space with the help of relation (10). As we mentioned, these integrals do not have a closed 
form solution, but fortunately they have an analytical solution in the limit of f3 — > (which is a good 
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approximation, because of the relation j3 



and we calculate them in terms of Fermi-Dirac and 



Bose-Einstein functions for general power law energy-momentum relation e = r]P a in Appendix B. 
The number of excited particles and pressure for fermions and bosons (see relations (B-5) and (B-8)) 
are given by 



N, 



^(T/ V f a T(3/a) h 3/a (z) 



1 - 3/3 (T/r?) 



a/- r(5/a) h 5/a {z) 
r(3/a) /i s/a (z) 



+ 6/3 2 (T/r?) 



4/a r(7/a) h 7/c » 



T(3/a) h {z 



+ ... 



(34) 



and 



V 
T 



^L(T/ V f a T(3/a) (z) 



2/ 



a T(5/a)/ i|+i (z) 
T(3/a)^ (*) 



+ 6/3 2 (T/n) 



4/c 



r(7/a)^ +i (z) 
r(3/a)^ +i (z) 



+ 



(35) 



where g is a weight factor showing the internal degree of freedom which for bosons and spin-1/2 
particles is equal 1 and 2 respectively. Here, h [/ (z) is given by relation (B-4) and reduces to Fermi- 
Dirac and Bose-Einstein functions according to the type of particles. Now, we can obtain energy of 
the systems by using of the definition (28) as 



I = ^(Th)' r(3/a) h, +i (z) 



a a 



h 3 



5f3(T/ v ) 



T(5/a)h §+i (z) 
T(3/a)hl (z) 



+ U[3 2 (T/rj) 



:a T(7/a)h i+i (z) 
T(3/a)h, I (z) 



+ 



(36) 



We can combine relations (35) and (36) to get equation of state for power law energy-momentum 
statistical systems 



V 

P 



T(5/q)/i 5 (z) 



T(3/a)h (z) 



■T{h/a)h %+i (z) 
T{3/a)h 3 (z) 



T(7/a)h l+i {z) 
T{3/a)h 3 (z) 



+ 



7T + 1 



(37) 



where p = U /V is the energy density. The l.h.s. of above equation denotes the equation of state 
(EoS) parameter uj which is a function of temperature in GUP framework. The first term in the 
r.h.s. is an ordinary EoS parameter and the other terms exhibit the quantum gravity effects which, 
in a cosmological point of view, are very small in the late times of cosmic evolution. However, in the 
planck scale they become important and play a determinant role in early universe thermodynamics. 
Then, we rewrite the above equation in the form 



u(/3,T) = oj 



r(5/a)h. (z) 



T(3/a)/i |+i (z) 
■ /r(5/a)/i |+i (*)y T a (?/ a )h l+i (z 



T(3/a)h (z) 

rv ' 



T{3/a)h (z) 



+ ... 



(38) 
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Figure 1: The effects of the minimal length on the EoS parameter to. As we see in the figure, the Eos parameter is 
a function of temperature in the GUP framework and get a negative value if temperature becomes more than critical 
temperature T c = {^)^T pl . 



Now, at the first, we consider above relation for non-relativistic particles. The energy-momentum 
relation for the non-relativistic particles is given by relation e = where m is the mass of the 
particles [a = 2 , 77 = ^— ) . In this case we have 



w(T) ~ 



2 r 



1 + 5m 



T 



(39) 



T pl , the second term in the r.h.s. of this relation 



It is clear that even at the Planck temperature T 
has a small value because of the large value of the Planck temperature T pl ~ 1.221 x 10 19 Gev. But 
for relativistic particles the energy-momentum relation is given by e = P (a = 1 = rj) and relation 
(38) for the bosonic gas in the limit z — > 1 becomes 



w(T) 



C(6) / T \2 
1 + 24-^4 ( * 



C(4) V Tj 



1440 



06) /C(6) _ 2 C(8)\1 (T_\ 
C(4) V C (4) C(6)^ \T P J 



(40) 



where (,{y) is the Riemmnn zeta function . Unlike the case of non-relativistic particles, the above 
relation does not converge in the limit of high temperature and all terms in the r.h.s. become 
important. Since the zeta function decreases with its argument, the coefficients of all the terms in 
r.h.s. expect the second one are negative. In this sense we offer the simplest temperature power law 
for the relation (40) as 



(41) 

where a and q > 2 are some constants which may be fixed by cosmological observation. An interesting 
feature of relation (41) is when the temperature increases and becomes larger than critical temperature 

T c = ia) 2q T Pl , for which we have a quintessence scenario for boson gaseous system. If we set T c = T pl 
then relation (41) describes an inflationary scenario, in figure 1 we have plotted the behavior of 00 for 
three values of q in which the dashed line shows the ordinary case and three curves show the effects 
of minimal length on cosmological parameter w as a function of temperature. It is clear that in the 
limit of low temperature all the curves tend to ordinary case which shows that the quantum gravity 
effects only become significant in high temperature i.e. in the Planck scales. 



1 - a 



T \2q 



T„ 
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4.1 Bose-Einstein condensation 



After set up the Bose-Einstein statistics, it is natural to consider famous phenomenon Bose-Einstein 
condensation. The equation (34) for bosons gives 



X = -^7 (T/rjf /a r(3/a) g (z) 



a 



h 3 



i - 3P(T/ V ) 



a/- E(5/a) g 5/Q (z) 
r(3/a) 5 3/a («) 



+ 6/3 2 (T/77) 



4/ 



a r(7/a) g 7/a (z 



r(3/a) <? 3/ J 



+ 



(42) 



We know from the ordinary case that the fugacity of the bosonic systems has a finite rang < z < 1 
and Bose-Einstein condensation occurs in the limit of z — > 1. Then, in this limit the above relation 
becomes 



N 



N 



i - 3p(T/ v ) 



i/q r(5/q)C(5/a) 
T(3/a)C(3/«) 



+ 6/3 2 (T/r?) 



4 /a r(7/a)C(7/a) 
r(3/a)C(3/a) 



+ ... 



(43) 



where iV 



= 4* CW° r(3/a) C(3/a) is the maximum number of excited particles in ordinary 
Bose-Einstein statistics. It is clear that the remaining terms in the r.h.s. are the quantum gravity 
effects and are small in low temperature. Therefore, for the non-relativistic particles with energy- 
momentum relation e = we get 



(44) 

Thus, the number of excited particles is smaller than the ordinary case and Bose-Einstein condensation 
occurs sooner in comparison with the ideal gas. However, the correction term due to the GUP effects 
is very small even in Planck Temperature. For relativistic particles with energy-momentum relation 
e = P, relation (43) becomes 



AT 



(2vrmT)i 
h3~ 



C(3/2) 



1 — 6m ■ 



K ^ F T 3 C(3) 



/ T \2 / T \4 

1-31 - +1812 - +.. 

^ pi pi 



(45) 



The above relation does not converge in high temperature and all the correction terms in the r.h.s. 
are important. 



4.2 Degenerate gas 

One of the interesting cases in Fermi-Dirac statistics is a completely degenerate gas which corresponds 
to the limit T — > 0. In this limit the mean occupation number < n e > in relation (16) takes a simple 
form as 



< n £ >= 




e > no, 



where /j,q denotes to the chemical potential in the limit T — > 0. Then, mean occupation number 
becomes independent of the temperature and acts as a step function. Now, the relations (17) and 
(18) become very simple and have an exact solution also in GUP framework. Then, relation (17) 
takes the form 
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N = - f a(e)de. 



(46) 



Here, there is a limited value for energy which is called Fermi energy e F , in terms of which one can 
also define the Fermi momentum P F . Upon using of relation (10), we have 



_ 4irgV f P F P 2 dP nVg 

~ ~h3~ Jo (1 + /3P 2 )3 ~ ~W 

which denotes to the number of particles in completely degenerate gaseous system with any energy- 
momentum relation in GUP framework. Up to the first order of B from the above relation one 
obtains 



arctan(V^-P F ) P F (PP F 2 - 1 



+ 



8(1 + BP, 



2\ 2 



(47) 



-(3(P F f-(P F f + P 0F 3 = 0, (48) 

where P 0p = ( ^y g ) 1 h is the Fermi momentum in ordinary case. We get P F = P 0p in the limit 
8 = 0, so we expand the above relation in the small but nonzero limit of GUP parameter 8 — > to 
get 



P ~ 

F 



3 o 24 , , 
1 + -/ P " ~ 175^ 



+ 



(49) 



The first term in r.h.s. of the above relation denotes the ordinary Fermi momentum and the second 
term is the most important effect of GUP on the Fermi momentum. Then, independent of that 
the gaseous system being relativistic or not, minimal length effects increase the value of the Fermi 
momentum and consequently increase the Fermi energy. 



5 Conclusions 

There are some evidences which support the idea of the existence of minimal length from quantum 
gravity theories. GUP theory imposes minimal length in quantum mechanics by modifying uncer- 
tainty relations, then the deformed Heisenberg algebra governs the phase space instead of the usual 
Heisenberg algebra. The customary procedure for exerting GUP conditions to a physical systems is 
finding the generalized coordinates and momenta which satisfy the deformed Heisenberg algebra. For 
the simplest choice, the Schrodinger equation becomes non-linear or it becomes a higher order differ- 
ential equation and it is not easy to solve it. In this paper we have extended statistical mechanics in 
GUP framework without fixing coordinates and momenta. We investigated the possibility of existence 
a general transformation in phase space which changes the usual Heisenberg algebra to a deformed 
one. Such a transformation changes the density of states in phase space. We obtained a deformed 
density of states for GUP conditions by the use of the Jacobian of the transformation. In this picture, 
effects of minimal length only changes the structure of phase space, so we conclude minimal length 
influences the situations of microstates in quantum scales and decreases the number of accessible 
microstates of the system. The advantage of this procedure is that we do not need to solve the mod- 
ified Schrodinger equation, instead we can work with the classical Hamiltonian together deformed 
density of state. In this method, we do not use the wave functions, so the exact solution of modified 
Schrodinger equations reserve their importance. Also we showed the non-commutativity of coordi- 
nates due to a GUP condition does not contribute in the density of states till the momenta commute. 
We develop the ensemble theory by this approach and we study the Maxwell-Boltzmann, Fermi-Dirac 
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and Bose-Einstein statistics. We found an interesting phenomenon in Maxwell-Boltzmann statistics: 
there is an upper bound for the number of excited particles for gaseous system (ideal gas and extreme 
relativistic gas) in the limit of high temperatures, which means that we have a condensation due to 
a GUP theory. We also found an interesting temperature dependent equation of state parameter oj 
for the bosonic gas based on which our analysis showed these effects are negligible for non-relativistic 
particles even in Planck scales. On the other hand, these effects become very important for the rela- 
tivistic bosonic gaseous systems and lead to an inflationary scenario in the limit of high temperatures. 
We showed there is a critical temperature which can be determined in the limit of the acceleration 
phase, but there is two free parameters that can be fixed with the cosmological observations. We 
studied two interesting full quantum cases, Bose-Einstein condensation and the completely degener- 
ate gas in the GUP framework. Our analysis shows the number of excited particles becomes smaller 
and the condensation occurs sooner. For the degenerate gas, we found a modified closed form for 
the number of particles in GUP framework and we obtained a modification to a Fermi energy and 
momentum. Finally we concluded all the corrections due to a minimal length are small in ordinary 
temperatures and will become important on the high temperature limit. 



A Jacobian in terms of Poisson brackets 

In this Appendix, we would like to show that it is possible to define general transformation in phase 
space which transforms coordinates Xi and momenta pi to generalized coordinates and momenta 
Pi as new phase space variable satisfy all the relations (5), (6) and (7). It is well-known from classical 
physics that the coordinates Xi and momenta pi satisfy the Poisson algebra {xi,pi} = 5ij. However, 
Darboux theorem implies that it is always possible to find variables Xi and Pi as a function of x% and 
Pi which satisfy relation 



{Xi,Pj} = f(Xi,Pj). (A-l) 

Therefore, the relation (6) is a particular case of the Darboux theorem and in GUP framework 
f(X,P) = f(P) = ih(l + (3P 2 ). On the other hand, relation (7) does not exert any constraint to 
our study. Thus, we can deduce that the relations (6) and (7) are valid simultaneously in general. 
But we want to have a transformation which satisfy also relation (5). We suppose general forms of 
non-commutativity as 



{X,Xj} =g{Xi,Pj), (A-2) 

where g is the arbitrary function of variables Xi and Pj. Again, it is possible to find variables Xi and 
Pj as a function of Xj and pi which satisfy relation (A-2), but we should have variables Xi and Pj as 
a function of Xi and pj which satisfy relations (5), (6) and (7) simultaneously. In general, we cannot 
prove that it is always possible to define general transformation which changes ordinary variables Xi 
and pj (which satisfy Poisson algebra) to GUP variables Xi and Pj (which satisfy deformed Poisson 
algebra). However, at least there is one choice Xi and Pi as 

Xi = Xi(l+ (3p 2 ) , 
Pi = Pi, 

which satisfy all the relations (5), (6) and (7) simultaneously. This implies to define transforma- 
tion between usual coordinates (xi,pi) and GUP variables pQ,Pj). The Jacobian of transformation 
can be expanded in terms of Poisson brackets in 2A^-dimensional phase space as [5] 



d{x t , p: 

d{xi , pi) 



2M 

E 



{Jin J12} ■■■ {J12M-11 ^ia/v}' 



(A-3) 
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where e denotes the Levi-Civita symbol and Jj represents the phase space variables as for odd i 
is a coordinate and for even i is a conjugate momentum Pj. In the above relation the Poisson 
brackets which denote only coordinates commutator (5), always multiply by Poisson brackets that 
purely denote the momentum commutator (7). Thus the coordinates Poisson brackets does not 
contribute because the momenta commute. Then only the Poisson brackets include both coordinates 
and momenta (6) contribute in relation (A-3) and Jacobian simplify to (9). 



B particle density and pressure in GUP framework 

The number of excited particles in GUP framework can be obtained from relation (19), by substituting 
mean occupation number < n £ > from (16), we get 



where a(e)de is the density of states in GUP framework determined by (10). Evaluation of this 
integral depends on the energy-momentum relation which unfortunately, there is no closed form 
solution for familiar power law energy-momentum. So we expand a{e)de for small f3. Up to the 
second order of /3 in three dimension (10) we obtain 

a(e) de ~ ( 1 - 3/3 P 2 + 6/3 2 P 4 + ...) . (B-2) 

We consider general power law energy-momentum relation e = r]P a which reduces to the ultra- 
relativistic case for rj = 1 = a and get the non-relativistic energy-momentum relation for r] = l/2m 
and a = 2. Then by substituting (B-2) in (B-l) results 



N 
V 



n 



4:71 (rrl \ 3,a 



a 



[<*> y a dy a/a [<*> y a 

Jo F^T^- 3/?(TA?) Jo Fv 



v ev + e 
+ 



y a ~ dy_ 
+ e 



y a ~ dy_ 
+ e 



+ 



(B-3) 



where n is the particle density and y = e/T. Integrals in the above relation are nothing but Fermi- 
dirac and Bose-Einstein functions for e = 1 and e = — 1 respectively. We define function h v (z) 
as 



h(z) = ^- r V dV = L{Z) 6 = +1 (B-4) 

where T(u) is a factorial function. Now particle density in terms of Fermi-Dirac and Bose-Einstein 
functions becomes 



^^(T/ V f a T(3/a) h 3/a (z) 



ah 3 



3/3 (T/ V ) 



2/ 



(x T(5/a) h . 



r(3/a) h 3/ Jz) 



+ 6/3 2 (T/n) 



Va T(7/a) h (z) 



r(3/a) h 3/a (z) 



+ ... 



(B-5) 



On the other hand, pressure of the system will be determined by relation (18). Again, for large volume 
we can replace summation by an integral. By using of relation (B-2) for general energy-momentum 
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relation we get 



V _ 4vr 

T eah 3 n 3/c 



/ ln(l +eze~ £ / T )e^ 1 de-3^ 2/a / ln(l + eze'^e"' 1 de 
Jo Jo 

+ 6/3 2 r? / ln(l + eze- £ / T )e a de + . 
Jo 



(B-6) 



Integrating by part and set y = e/T yields 



V 4vr 3 /a 



00 y 3/a dy 
z-V + e 

18 



x y° /a dy 



e" + e 



18 2 4/a r 00 y^dy 
7 Jo ^ e + e 



(B-7) 



The above integrals are again Fermi-Dirac and Bose-Einstein functions and we can express them in 
terms of h v (z) as 



f = %> (T/7])3/a r(3/a) 



1 - 3/3(T/77) 



2/0 



r(5/a)/ l|+i (z) 
T(3/a)h"(z) 



+ 6/3 2 (r/77) 



ia r(7/a)^ +i (z) 
r(3/a)/i a+i (z) 



7T + 1 



+ 



(B- 



Relations (B-5) and (B-8) are essential for our study in thermodynamics of bosons and fermions in 
GUP framework. 
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